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Abstract. The giant radio galaxy M87 is usually classified as a Fanaroff-Riley class I source, suggesting that M87 is 
a mis-aligned BL Lac object. Its unresolved nuclear region emits strong non-thermal emission from radio to X-rays 
which has been interpreted as synchrotron radiation. In an earlier paper we predicted M87 as a source of detectable 
gamma ray emission in the context of the hadronic Synchrotron-Proton Blazar (SPB) model. The subsequent 
tentative detection of TeV energy photons by the HEGRA-telescope array would, if confirmed, make it the first 
radio galaxy to be detected at TeV-energies. We discuss the emission from the unresolved nuclear region of M87 in 
the context of the SPB model, and give examples of possible model representations of its non-simultaneous spectral 
energy distribution. The low-energy component can be explained as synchrotron radiation by a primary relativistic 
electron population that is injected together with energetic protons into a highly magnetized emission region. We 
find that the 7-ray power output is dominated either by /i^/vr* synchrotron or proton synchrotron radiation 
depending on whether the primary electron synchrotron component peaks at low or high energies, respectively. 
The predicted 7-ray luminosity peaks at ~100 GeV at a level comparable to that of the low-energy hump, and this 
makes M87 a promising candidate source for the newly-commissioned high-sensitivity low-threshold Cherenkov 
telescopes H.E.S.S., VERITAS, MAGIC and CANGAROO III. Because of its proximity, the high-energy spectrum 
of M87 is unaffected by absorption in the cosmic infrared (IR) background radiation field, and could therefore 
serve as a template spectrum for the corresponding class of blazar if corrected for mis-alignment effects. This 
could significantly push efforts to constrain the cosmic IR radiation field through observation of more distant 
TeV-blazars, and could have a strong impact on blazar emission models. If M87 is a mis-aligned BL-Lac object 
and produces TeV-photons as recently detected by the HEGRA-array, in the context of the SPB model it must 
also be an efficient proton accelerator. 

Key words. Galaxies: active - Galaxies: individual: M87 - Gamma rays: theory - Radiation mechanisms: non- 
thermal 



1. Introduction 

The FanarofF-Riley (FR) class I giant radio galaxy M87, situated nearly at the center of the Virgo cluster, was the 
first extragalactic jet to be discovered (Curtis (1918)), and has since then been intensively observed at all wavelengths. 
Its proximity (~ 16.3 Mpc; Cohen et al. (2000)) makes it an interesting laboratory for testing and understanding 
extragalactic jets of radio-loud Active Galactic Nuclei (AGN) and their powering engines. 

Because M87 is sufficiently near for ultra-high-energy cosmic rays (UHECRs) to be little affected by the GZK- 
cutoff at ~ 5 X 10^^ eV (Greisen (1966), Zatsepin & Kuz'min (1966)), and because its size scales could allow magnetic 
confinement of the most energetic cosmic rays (Hillas (1984)), M87 has long been considered as one of the prime 
candidate sources of high energy cosmic rays. This idea has recently received some support by the suggestion that 
possible clustering observed in the arrival directions of the UHECRs can be understood in terms of deflection of 
UHECRs from M87 by our Galaxy's magnetized wind assuming a Parker-spiral magnetic structure (Ahn et al. (1999); 
Biermann et al. (2001)). It appears, however, that such deflection in the Galactic wind may be insensitive to the 
direction of the cosmic ray sources (Billoir & Letessier-Selvon (2000)). Nevertheless, Protheroe et al. (2003) found 
by using the cosmic ray output predicted in the SPB model that M87 could explain the observed UHECR flux if 
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the magnetic field topology between M87 and our Galaxy were favourable. They found that UHECR with energies 
above 10^° eV could easily be produced because neutrons produced in the pion photoproduction process would be 
relativistically beamed along the jet direction and Dopplcr boosted in energy. Even though M87's jet is mis-aligned 
with respect to our line-of-sight, these Doppler boosted neutrons escape from the jet and decay into UHECRs which 
maintain their Doppler boosted energies and may propagate in all directions, including towards our Galaxy if the 
magnetic field topology were favourable. Of course we note that M87's nuclear region is not the only possible source 
of the UHECRs observed at Earth; see e.g. Protheroe & Clay (2004) for a recent review. 

According to the unification model of AGN (e.g. Urry & Padovani (1995)) FR I radio galaxies, with their jet axis 
at a large angle to our linc-of-sight, arc the parent population of BL Lac objects whose jets arc closely aligned to 
our line-of-sight. This motivates us to consider M87's nuclear region as a mis-aligned blazar of BL Lac type. The 
spectral energy distribution (SED) of BL Lac objects can usually be explained satisfactorily by either leptonic or 
hadronic blazar emission models. Bai & Lee (2001) have discussed M87 on the basis of the leptonic Synchrotron-Self 
Compton (SSC) model where synchrotron photons produced by interactions of relativistic electrons with the ambient 
magnetic field serve as the target photons for inverse Compton scattering by the same electrons. By interpreting the 
non-thermal radiation from the radio to the X-ray band as synchrotron emission with luminosity peaking in the far- 
ultraviolet, the authors considered M87 to be a mis-aligned high-frequency peaked BL Lac (HBL), and predicted 7-ray 
emission with an inverse Compton peak at ~ 100 GcV. The predicted inverse Compton flux is consistent with the recent 
HEGRA detection of M87 (Aharonian et al. (2003), see Sect. 2). Detectable TcV-emission from Comptonization of 
galactic photon fields has recently been suggested by Stawarz et al. (2003). In contrast to former models, they consider, 
however, the large scale jet to be the site of 7-ray production. 

While in leptonic models a relativistic electron-positron plasma is usually assumed to be responsible for the non- 
thermal jet radiation, in hadronic models a relativistic proton-electron (p e~) plasma is assumed to be the main 
constituent of the jet material. In the hadronic Synchrotron-Proton Blazar (SPB) model, proposed recently by, e.g., 
Miickc & Protheroe (2001), accelerated protons interact with the synchrotron radiation field produced by the co- 
accelerated electrons via meson photoproduction and Bethe-Heitler pair production and, more importantly, with the 
strong ambient magnetic field emitting synchrotron radiation (mesons and muons also emit synchrotron radiation). 
The SPB model neglects external photon field components, and this seems appropriate for BL Lac objects and their 
parent population which possess only weak accretion disks. Miicke & Protheroe (2001) have shown that this model 
can reproduce the commonly observed double-humped blazar SED. Hadronic models require high proton energies that 
can only be achieved in a highly magnetized environment where synchrotron losses can become severe. Magnetic field 
values around lO'^ G are thought to exist near the horizon of a supermassive black hole (Blandford & Znajek (1977)) 
with a mass of -^10^ Mq as estimated for M87 (Marconi et al. (1997)). However, with M87's rather low accretion 
rate if the equipartition value of B scales with M, and assuming magnetic energy flux conservation, magnetic fleld 
strengths of order 10-100 Gauss are expected within 30 Schwarzschild radii rg where the jet is probably formed 
(Junor & Biretta (1999)). In the present work, we discuss in more detail than our earlier work (Protheroe et al. (2003)), 
and in the context of the recent HEGRA detection (Aharonian et al. (2003)), the nuclear (core) emission, i.e. from 
the M87 jet, in the framework of the SPB-model. 

In Section 2 we summarize the data on M87's core emission. In section 3 we give a brief model description, and 
calculate the steady-state synchrotron component as described in the appendix. The modeling procedure is described 
in section 4, and we conclude with a summary and discussions in section 5. 

2. The data 

Speculation that M87 could be a powerful accelerator of cosmic rays triggered space-based 7-ray detectors and ground- 
based high-sensitivity Cherenkov telescopes to search for 7-ray emission from this radio galaxy. Until recently, only 
upper limits were obtained. From EGRET data Reimer et al. (2003) obtained F(>100 MeV) < 2.2x10"^ cm"^ 
s-^ Using Whipple data from 2000-2001 (Le Bohec et al. (2001)) and 2002-2003, Le Bohec et al. (2003) obtained 
F(>250 GeV) < 2.6 x 10~" cm'^ s'^ Using 1998-1999 data from HEGRA telescope array Getting et al. (2001) 
obtained F(> 720 GeV) < 1.45x10"^^ cm~^ s~^ (at the 2>a level). By doubling the data set and applying a more 
sensitive analysis method, the HEGRA team has recently been published the first (though tentative) detection of > 730 
GeV photons at the 4(7 level (Aharonian et al. (2003)). This detection, which does not contradict the Whipple upper 
limits, places the first data point in the so-far rather unconstrained high energy regime of M87's SED, and has motivated 
us to refine our previous modeling attempt and predictions for 7-rays from this source (Protheroe et al. (2003)). 

The non-simultaneous SED is shown in Fig. 1 and 3. The data imply that the synchrotron spectrum from the 
primary electrons must exhibit a break frequency at either a few lO^^Hz (Perlman et al. (2001)) or around lO^^Hz. 
Break frequencies above lO'^'^Hz as proposed by Bai & Lee (2001) seem unreasonable to us given the fact that, firstly, 
the EUVE data point from Berghofer et al. (2000) with its considerably lower resolution compared to the other data 
points should effectively be considered rather as an upper limit, and secondly, the recent flux measurements from the 
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Chandra observatory (Wilson & Yang (2002)) point to a much steeper spectrum in the X-ray band than anticipated 
by Bai & Lee (2001). Note, however, that variability effects may play a crucial role here. Another interesting spectral 
feature is the strong steepening by Aa w 1 that occurs around IR/optical wavelengths, which can not easily be 
explained by a transition from escape dominated to synchrotron cooling dominated electron energy losses. 

Images from the VLA (Biretta et al. (1995)), HST (Sparks et al. (1996), Biretta et al. (1999)), Gemini 
(Perlman et al. (2001)) and Chandra (Wilson & Yang (2002). Harris et al. (2003)) of the so- far unresolved core of 
M87 are at present at the sub-arcsec scale, giving a metric resolution of order 1-10 pc. Sub-mas scale cm and mm 
wavelength intercontinental VLBI (Junor & Biretta (1995), Junor & Biretta (1999)) provide the highest linear reso- 
lution of ^0.01 pc achieved on any cxtragalactic radio jet so far, and thereby place the most stringent upper limit on 
the size of the radio emitting region. Striking variability in the core region has not only been observed in the radio to 
optical band, but has also been deduced from Chandra X-ray monitoring in 2002 (Harris et al. (2003)). An observed 
flux increase of about 20% which has been measured within 46 days can be transformed into a doubling time of about 
77 days, and provides a limit for the source size of the X-ray emitting region of i? ~ O.lpc. 

HST data (Biretta et al. (1999)) show that features within the first arcsec of the jet move only at sub-luminal 
speeds, while at larger distances from the core super-luminal motion is observed. In this work we assume that the 
pattern speed equals the flow velocity of the knots, or "plasmoids". The upstream knot closest to the core has an 
apparent speed of 0.63±0.23c which we use here to constrain the beaming factor for the nuclear emission. For a jet angle 
between 10° — 40° as suggested from VLA and HST proper motion studies (Biretta et al. (1995), Biretta et al. (1999)) 
we find that Doppler factors in the range D = 1.5-3 are consistent with the apparent bulk speed. 

3. The model 

We assume the emission region, or "blob" , in an AGN jet moves relativistically with Lorentz factor Tj and velocity 
/3jC along the jet axis. We further assume that rclativistic (accelerated) protons, whose particle density N'^ follows 
a power-law spectrum oc 7^ "" in the range 2 < 7^ < 7p^„iax (primed quantities are in the jet frame), are injected 
instantaneously into a highly magnetized environment {B' is constant within the emission region), and that they 
remain quasi- isotropic in the jet frame due to pitch-angle scattering. The proton energy-loss processes considered in 
the model are photomeson production, Bethe-Heitler pair production, proton synchrotron radiation and adiabatic 
losses due to jet expansion. Synchrotron radiation prior to their decay from tt^ (from photomeson production) and /i^ 
(from TT^ decay) becomes important in highly magnetized environments (Rachen & Meszaros (1998)), and is taken 
into account in our calculations. We assume that the maximum particle energies are limited by the balance between 
energy gain and loss rates. The acceleration rate for any acceleration mechanism is dE' /dt' = S,(E')ZecB' where 
< 1 is the acccilcration rate factor and Ze is the charge. If particles gain energy by diffusive shock acceleration 
and the spectra of both electrons and protons are cut off by synchrotron losses at Lorentz factors 7e,max and 7p,max, 
respectively, then 



if the diffusion coefficient has energy dependence k(E) oc E^ (sec Appendix A). Thus for a Kolmogorov spectrum of 
turbulence, for which (5=1/3, ^(£^e,max)/'C(£^p,max) = {rn,./mp)^ « 3 x 10"^. 

The rclativistic primary e~ , injected into the emission region with a power law particle distribution cx E~""' , radiate 
synchrotron photons that manifest themselves in the blazar SED as the synchrotron hump, and serve as the target 
radiation fleld for proton-photon interactions, and for the subsequent pair-synchrotron cascade which develops as a 
result of photon-photon pair production in the magnetized blob. The steady-state primary electron spectrum in the co- 
moving frame of the emission region is calculated as in Appendix A taking into account synchrotron and escape losses. 
The resulting synchrotron radiation from this particle distribution is then corrected for synchrotron-self absorption. 
For B' > 0.6(UpijQt/10^°eV cm~^)^/^ Gauss the target photon density Uphot smaller than the magnetic field energy 
density. Thus Inverse Compton losses can in most cases be neglected in the SPB model, and the corresponding SSC 
component is expected to be low compared to the primary synchrotron 'hump'. However, for completeness we have 
calculated explicitely the SSC radiation from the primary electron component (see Appendix A). 

The pair-synchrotron cascade redistributes the photon power to lower energies where the photons eventually escape 
from the emission region of size R' . The cascades can be initiated by photons from 7r°-decay ("tt" cascade"), electrons 
from the tt^ ^ /x^ decay ("tt''^ cascade"), p-synchrotron photons ( "p-synchrotron cascade"), charged /x-, tt- and 

/C-synchrotron photons ( "/x^-synchrotron cascade") and from the proton-photon Bethe-Heitler pair production 
("Bethe-Heitler cascade"). Direct proton and muon synchrotron radiation is mainly responsible for the high energy 
hump whereas the low energy hump is dominated by synchrotron radiation from the primary e~, with a contribution 
of synchrotron radiation from secondary electrons produced by the p-synchrotron and /i^-synchrotron cascades. The 
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contribution from Bethe-Heitler pair production turned out to be negligible. For our calculations we use a Monte-Carlo 
method and utilize the recently developed SOPHIA code for the photohadronic event generation (Miicke et al. (2000)). 
In practice, to save CPU-time the target photon field is parametrized as a multiple broken power-law, which is then 
used as an input into the SOPHIA code. 

4. Modeling the SED of M87 

From variability and direct imaging arguments an upper limit for the size of M87's emission region of a few lO^^cm 
in the observer frame can bo deduced. Together with the constraints from proper motion measurements, a reasonable 
parameter space for the modeling procedure is R' = 10^'' '^^'cm and bulk Doppler factors D = 1.5 . . .3. Furthermore 
we demand approximate equipartition between magnetic u'^ and particle energy density {u'^+u'^) « u'^. This in general 
minimizes the total jet power for a given parameter set (Miicke & Protheroe (2001)). The observed synchrotron hump 
in the SED implies a break in the primary electron spectrum. A peak in the "synchrotron hump" in the SED is 
expected at either mm-wavelengths or in the optical. In the following we shall model both possibilities in turn, i.e. a 
high and low energy peak. 

4.1. High-energy peaked synchrotron component (Model H) 

The HEGRA-detection at sub-TcV energies places an important constraint on the models: protons must be accelerated 
to energies above 10^°GeV. This can only be achieved by high magnetic field values, and/or a thin target photon field 
to prevent excessive losses at the highest energies. Typical magnetic field values in the SPB-model lie around several 
10 G. Fig. 1 shows a reasonable representation of the data where we have used _B' = 30 G (u^=2xl0^'^ eV cm~^) 
and the primary synchrotron component peaking at about 1 eV. With a Doppler factor of D = 2 and an observer- 
frame size of the emission region of ~ 10^^ cm the target photon density is Mphot ~ 3x10^'^ eV cm~^, and so the 
SSC component is negligible in Model H. Proton synchrotron losses dominate at the highest proton energies for this 
parameter set (see Fig. 2), and determine the cutoff energy and the 7-ray output from ~ 1 TeV down to ~ 10 MeV. At 
0.1-10 MeV synchrotron radiation from a secondary e^ population produced by the reprocessed /x^/tt^ synchrotron 
radiation dominates and produces a broad "valley" in the SED between the low and high energy humps that possess 
approximately equal power. An acceleration rate factor at the maximum proton energy ^(£^p^i„ax) ~ 1 is necessary to 
allow the injection proton spectrum to extend up to 7p^rnax = 3 x 10^°. 

The primary electron synchrotron spectrum shows a low-energy break at the synchrotron self- absorption turnover 
energy of ~ a few 10^'' eV, followed by a nearly flat power distribution, and then a turnover at about a few eV with 
a subsequent steep tail due to the cutoff in the electron distribution. Note that interpreting the strong steepening at 
a few eV in the data (sec Sect. 2) as being due to a cutoff in the electron spectrum allows spectral breaks larger than 
0.5. The high magnetic field leads to a dominance of synchrotron losses throughout the emitted low-energy component 
(the escape loss dominated energy range lies below the synchrotron-self absorption turnover frequency). 

The total jet power for this parameter choice is Ljot ~ 3xl0'*'^erg/s, about the value for the jet kinetic power 
derived by Reynolds et al. (1996), but still below the nuclear power of M87 for accretion at the Bondi spherical rate 
(~ 5 X 10*^ erg/s, Di Matteo et al. (2003)). The small radiative efficiency of the accretion disk and the fact that the 
radio power of M87 is low provides a natural explanation that M87 is a FRI source. Owen et al. (2000) have calculated 
the total bolometric luminosity of the order of 10^^ erg/s, in agreement with the total radiative output in the present 
models, which suggests that the jet in M87 is also a low efficiency radiator and that M87 is currently in a dormant 
activity stage. 

4.2. Low-energy peaked synchrotron component (Model L) 

Fig. 3 shows an example for a model in which the primary electron synchrotron component peaks at ~ 10~^eV. Here 
we have chosen the size of the emission region to be of the order of the limiting size from direct VLBI imaging, i.e. 
~0.01 pc (see Sect. 2). With D = 1.5 and R ~ 10^^ cm the jet-frame target photon density is low, Uphot ~ 2x10^ eV 
cm~^. A relatively low magnetic field strength of i?' = 5 G gives equipartition between magnetic and particle energy 
densities. 

Despite the lower target photon energy density (down by a factor ~10 compared to model H) pion photoproduction 
losses dominate at the highest energies for this parameter set (see Fig. 2) since the magnetic energy density is reduced 
by even more (down by a factor 36 compared to model H). In fact, even for models with identical magnetic and target 
photon energy densities, a lower break-energy in the target spectrum would cause a turnover in the n production 
losses at correspondingly higher proton energies and this would in turn result in a significantly higher tt production 
loss rate at the highest proton energies. This is demonstrated in Fig. 2 where we compare the tt production loss rate 
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Models (H) and (L) for 7r-photoproduction (tt, 
solid lines) and proton synchrotron radiation (p 
syn, dashed line). While in model (H) p syn- 
chrotron losses dominate, in model (L) tt produc- 
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distribution as in model (L). Lowering the target 
field's peak energy leads to an increase of the n 
production rate at the highest proton energies. 



of the present parameter set (Model L) with the corresponding rate for the same target spectrum but normalized up 
to the photon energy density as used in Model H. We conclude that in radiation fields that peak at high energies 
p-synchrotron losses dominate, while in low-energy-peaked target photon fields tt production losses dominate at the 
highest proton c^ucrgies for equal magnetic and target photon energy densities. 

For modeling the primary synchrotron spectrum peaking at around 10~^ eV we have injected a softer electron 
distribution {ue = 2.1) into the emission region. The turnover at this energy is due to synchrotron-self absorption 
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becoming dominant at radio wavelengths. A gradual steepening followed by a steep decline occurs above ~ 1 eV caused 
by the cutoff in the electron distribution. Again synchrotron losses dominate the steady-state electron spectrum above 
the self absorption turnover energy. However, because the ratio Uphot/'^B significantly higher than in Model H, the 
SSC component (solid curve at X-ray energies in Fig. 3) is potentially important in Model L. 

Because pion photoproduction losses dominate over proton synchrotron losses for the present parameter set, the 
predicted 7-ray spectrum above 10 GcV is determined by synchrotron emission by charged pions and muons rather 
than by the protons themselves. Muon/pion synchrotron radiation naturally produces a high energy hump extending 
even up to TeV energies. Because pions and muons have a lower rest mass than protons, their synchrotron emission 
peaks at higher photon energies than the proton synchrotron radiation for the same Lorentz factor and magnetic 
field. While /U^/tt^ synchrotron radiation dominates the energy output around the peak energy at ~100 GeV, the 
emission in the EGRET energy range is due to proton synchrotron radiation, and extends down to soft 7-rays where 
the synchrotron radiation from the 7r-cascades takes over. Again, acceleration rate factor at the maximum proton 
energy ^(i^p^max) ~ 1 is needed to explain photon emission up to TeV-energies. 

Models with even lower jet-frame target photon densities could also be consistent with the data (except models with 
very low magnetic field values < 2 G). However, the energy density stored in particles would be orders of magnitude 
below the magnetic energy density. If the size of the emission region were of order lO^^cm or less, the resulting target 
photon density would reach values above 10^^'^eV/cm^, and photopion production losses would cut off the injected 
proton spectrum at ^ lO'^GcV. In such models TcV-emission at a flux level as detected by HEGRA would be difficult 
to explain by proton or fi^/w^ synchrotron radiation. 

5. Summary and discussion 

We have made SPB model flts to the non-simultaneous SED of M87's nuclear emission, and find that all parameter sets 
which satisfactorily represent the data predict the main contribution to the high energy luminosity at about 100 GeV 
to be due to either /i^/7r^ synchrotron or proton synchrotron radiation depending on whether the primary electron 
synchrotron component peaks at low (Model L) or high (Model H) energies, respectively. 

In the EGRET energy range, the lower synchrotron peak energy model (Model L) predicts a softer spectrum than 
the fit with a higher synchrotron peak energy (Model H). While it is obvious that EGRET's sensitivity was more 
than an order of magnitude above the expected flux level from M87, we find that the satellite-based 7-ray instrument 
GLAST might possibly detect a weak signal from this radio galaxy (see Fig. 4). 

In all the fits we have presented, the high energy radiation cuts off with a strong steepening in the TeV range in 
agreement with the spectral limits from the HEGRA observation. We also find that the HEGRA detection at >730 
GeV can only be explained if the proton acceleration rate is extremely high (C(-E'p,max) ~ !)• We therefore expect M87, 
if it is indeed a mis-aligned SPB, could be an important source of UHECRs (see also Protheroe et al. (2003)). 

For almost all proposed models of particle acceleration in different astrophysical environments, £,{E) remains a 
rather uncertain model parameter. On the other hand, any postulation of acceleration of high energy protons in 
compact 7-ray production regions actually implies that ^(-Ep^inax) these energies should be close to unity, which 
corresponds to the maximum theoretically possible acceleration rate based on simple geometrical consideration (e.g. 
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dashed-triple dotted lines; H.E.S.S. phase I: solid lines; MAGIC (from: "The MAGIC Telescope Project Technology and per- 
formance aspects": http://hegral.mppmu.mpg.de/MAGICWeb): dashed-dotted lines; CANGAROO III (from: "Status Report 
of the CANGAROO-III Project": http://icrhp9.icrr.u-tokyo.ac.jp/c-iii.html): dotted lines. For the GLAST sensitivity (from: 
http://glast.gfsc.nasa.gov/resources/brochures/gsd/) (long dashed line) a source photon spectrum oc E~^'^ is used. 



Hillas (1984)). An interesting possibility could be particle acceleration at the annihilation of magnetic fields in the fronts 
of poynting flux dominated jets (Blandford (1976); Lovelace (1976)). It has been argued that this mechanism could 
provide effective acceleration of extremely energetic protons with ^(i?p „,ax) ~ 1 (Haswell et al. 1992). A quantitative 
investigation of particle acceleration mechanisms is beyong the scope of this paper. 

Both model fits presented seem to under-predict the emission in the radio domain as compared to the observations. 
Our modeling, however, assumes the same size for the emission region at all energies, while the data indicate a smaller 
width of the optical than the radio jet (Sparks et al. (1996)) though with roughly the same morphology. In addition, 
the inter-knot region is observed to be weaker in the optical than in the radio band (Sparks et al. (1996)). It appears 
therefore reasonable to attribute the missing flux in our model to the inter-knot region, and to a larger blob size in 
the radio band as compared to higher frequencies. 

As previously noted, the Chandra data lie above the extrapolation of the optical spectrum to higher energies 
(Wilson & Yang (2002)). We point out that the modelled SED of M87 is based on non-simultaneous data, and that 
the X-ray flux could have been much lower than the Chandra data suggest at the time of the optical observations. 
Another critical point is that different resolutions of the images at the various energies have been used in the literature 
for the flux determinations. We have used for our compilation of the SED the highest resolution data available at each 
frequency, ranging from arcmin (HEGRA) to sub-arcsecond (Chandra, Hubble, VLBI) scales, and could introduce 
additional non-negligible uncertainties into the flux measurements. However, a flatter X-ray spectrum consistent with 
the Chandra data might be achieved if either the magnetic fleld is highly inhomogenous, or a secondary e^ population 
is responsible for the X-ray flux (Wilson & Yang (2002)). 
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While, in principle, the SPB model provides secondary synchrotron emitting e s from the various cascades which 
may extend even into the X-ray domain, the fits presented can not easily explain the high flux level observed at these 
energies as electron synchrotron radiation or as part of the high energy hump cascade component. Nevertheless, for 
magnetic fields of order a few Gauss SSC radiation might become detectable also in hadronic models. Fig. 3 shows 
that the SSC component peaks at X-ray energies with a spectral signature that is in agreement with the Chandra 
observations, but its flux is roughly an order of magnitude too low to explain the Chandra data. If the magnetic field 
or the size of the emission region were a factor ~ 3 lower model L would predict the Chandra X-rays as SSC photons. 
The latter change moves model L slightly in the direction of model H. Thus, X-ray variability might be related to 
relatively small changes in the model paramaters, raising or lowering the importance of inverse Compton losses with 
respect to synchrotron losses. 

Another possible source for the observed emission in the Chandra band might be a contribution, either directly 
or reprocessed through cascading, from photon fields other than the primary electrons' synchrotron radiation. The 
spectral continuum data do not show any thermal bump from the putative accretion disk. Di Matteo et al. (2003) have 
shown that an advection dominated accretion disk could account for a large fraction of the observed X-ray nuclear 
flux. The radiative efficiency is extremely low, but the accretion rate is found to be large enough for Comptonization 
of the synchrotron emission of the disk or the thermal bremsstrahlung emission to dominate the X-ray emission 
(Di Matteo et al. (2003)). A sudden drop of one order of magnitude of the accretion rate could lower the X-ray disk 
output by --^ 2 orders of magnitude. 

For an advection dominated disk with a X-ray luminosity ^disk ~ 7 X 10^° erg/s (Di Matteo et al. (2003)) the total 
energy density in the jet frame Mdisk be derived through the transformation (see e.g. Dermer & Schlickeiser (2002)) 

where r is the distance of the jet plasma blob from the black hole, /i(r) = r/(r^ -|- rf^Y^'^ and we have approximated 
the disk's radiation field as a luminous ring of radius rjn that illuminates the moving blob. Assuming a inner accretion 
disk radius of rin = 1.23rg (for an extreme Kerr black hole of mass M = lO^M©) and taking Tj = 1.5 one obtains 
"disk 2 • 10^°eV cm~^ at r = lO^^cm and ~ 2 • lO^eV cm~^ at r = lO^^cm. Hence, for r lO^^cm the accretion disk 
radiation proves to be unimportant as a target field for cascading and photon-particle interactions in M87 compared 
to the primary electron synchrotron emission. 

Recently Donea & Protheroe (2003) have constrained the torus temperature of the torus to < 250 K using existing 
data from the literature. On the other hand, during an extreme flaring state, related to the accretion rate changes or 
to a spin flip of the central black hole, the torus could undergo enough heating to become 'visible'. This alters the high 
energy part of the spectrum above several hundreds of GeV, as is discussed in Donea & Protheroe (2003). Therefore, 
a visibility-state of the torus (if present) could be achieved at the cost of not being able to observe very high energy 
gamma rays from the nucleus of M87. Regular monitoring of M87 at VHE gamma-rays and IR frequencies could be 
important to elucidate the problem of existence or non-existence of a dusty torus in M87. For a temperature of the 
torus radiation of < 250K the co-moving frame energy density is 2 • lO^F^ cV cm~^ <C Wphot M87, and is therefore 
negligible as target photon field. The star and dust contribution of M87's host galaxy has been estimated to 630F^ eV 
cm~^ and 6.3r| eV cm~^ in the jet frame, respectively (Stawarz et al. (2003)), and can obviously also be neglected 
regarding M87's synchrotron radiation density of order 10^° eV cm~^. 

Fig. 4 shows that the recently commissioned Cherenkov telescope array VERITAS, the southern arrays H.E.S.S. and 
CANGAROO III (though at large zenith angles > 45°), and MAGIC may be able to detect M87. The predicted integral 
fluxes > 100 GcV for Models H and L arc ^ 4 x 10"^^ cm~^ s~^ and ~ 4 x 10~^^ cm~^ s~^, respectively. We have used 
A. Konopelko's simulator for the H.E.S.S. response (http://pluto.mpi-hd.mpg.de/ konopelk/WEB/simulator.html) to 
estimate the necessary observation time for statistically signiflcant detection. A 10 h on-source observation with the 
full phase I (four telescopes) H.E.S.S. array would result in a 8 — 9(J detection (expected cosmic ray rate is ~ 0.7 s"'^, 
7-ray rate is 0.055 s~^) in the case of a high-energy peaked photon target, and 4 — 5cr detection in the case of the 
low-energy peaked photon target for 300 h of usable data assuming the source at zenith (expected cosmic ray rate is 
here ~ 0.7 s~\ 7-ray rate is 0.006 s~^). Since the sensitivity of VERITAS (Weekes et al. (2002)) is similar to that 
of H.E.S.S., similar numbers can be expected for VERITAS observations. In Fig. 4 we summarize the minimum flux 
for a 50 h observation (with statistics exceeding 10 photons and a signal detection at a level of at least 5(t) using 
the phase I H.E.S.S. array, the VERITAS array, CANGAROO III and MAGIC (assuming the source at zenith) and 
GLAST, in comparison the the predicted high energy fluxes. Note, however, that these predictions are based on a 
non-simultaneous observed SED and, depending on the actual activity state of M87, the predicted fluxes and spectra 
may change signiflcantly. In addition, absorption of 7-rays in infrared radiation from a putative torus could affect the 
spectrum above 1 TeV if the torus temperature Ttorus were higher than 250 K, and above 200 GeV if Ttorus > 1000 K 
(Donea & Protheroe (2003)). Work is in progress to make SPB model fits to other nearby FR I radio galaxies. 
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In both models presented here, the power output in the high energy hump is roughly equal to the power output 
in the low energy hump of the SED. Because of M87's proximity, absorption of sub-GeV/TeV-photons in the cosmic 
infrared background radiation field is not expected to affect the spcctnim below ^ 50 TcV. The observed spectral 
behaviour at high energies should he intrinsic to the source. Tracing the spectrum at GeV-TeV-energies would give a 
7-ray spectrum that for the first time includes an unabsorbed (by radiation fields external to the source) cutoff. These 
data could serve as a typical template BL Lac spectrum at source after correcting for M87's jet mis- alignment. By 
comparing this template with BL Lac spectra at high redshifts, meaningful constraints for the extragalactic background 
radiation field around IR wavelengths can be derived. 

Acknowledgements. We thank A. Konopelko for providing us with his H.E.S.S. response simulator and fruitful discussions. AR's 
research is funded by DESY-HS, project 05CH1PCA/6, and that of RJP and ACD by an ARC Discovery Project grant. 

References 

Aharonian, F. and the HEGRA Collaboration 2003, A&A 403, LI 

Ahn, E.-J., Medina- Tanco, C, Biermann, P.L., Stanev, T. 2000, astro-ph/9911123 

Bai, J.M. & Lee, M.G. 2001, ApJ, 549, L173 

Berghofer, T.W., Bowyer, S., Korpela, E. 2000, ApJ, 535, 615 

Biermann, P.L. & Strittmatter, P.A. 1987, ApJ, 322, 643 

Biermann, P.L. 1997, J. Phys. G, 23, 1 

Biermann, P.L., Ahn, E.-J., Kronberg, P.P., Medina- Tanco, C, Stanev, T. 2001, in Physics and Astrophysics of Ultra-High- 
Energy Cosmic Rays, vol. 576, ed. M. Lemoine & G. Sigl (Lecture Notes in Physics) 181 
Billoir, P., Letessier-Selvon, A. 2000, astro-ph/0001427 
Biretta, J.A., Sparks, W. B., Macchetto, F. 1999, ApJ, 520, 621 
Biretta, J.A., Zhou, F., Owen, F.N. 1995, ApJ 447, 582 
Biretta, J.A., Stern, CP., Harris, D.E. 1991, AJ, 101, 1632 
Blandford, R.D. 1976, MNRAS 176, 465 
Blandford, R.D. & Znajek, R.L. 1977, MNRAS, 179, 433 
Cohen, J.C. 2000, AJ, 119, 162 
Curtis, H.D. 1918, Pub. Lick. Obs. 13, 31 
Dermer, CD. & Schhckeiscr, R. 2002, ApJ, 575, 667 

Di Matteo, T., Allen, S.W., Fabian, A.C., Wilson, A.S., Young, A.J. 2003, AJ, 582, 133 
Donea, A.-C. & Protheroe, R.J. 2003, Prog. Theor. Phys. SuppL, 151, 186 

Drury, L.O'C 1983, Rep. Prog. Phys., 46, 973 
Falckc, H. & Biermann, P.L. 1995, A&A, 293, 665 

"GLAST: Exploring nature's highest energy processes with the Gamma Ray Large Area Space Telescope" NASA document 
NP-2000-9-107-GSFC, February 2001, p29. http://glast.gfsc.nasa.gov/resources/brochures/gsd/ 

Getting, N. and the HEGRA Collaboration 2001, in Proc. 27th Int. Cosmic Ray. Conf., Hamburg, vol. 7, ed. M. Simon et al. 
(Copernicus GcscUschaft, Katlcnburg-Lindau) 2669 

Greisen, K. 1966, Phys. Rev. Lett., 16, 748 

Hillas, A.M. 1984, ARA&A, 22, 425 

Harris, D.E., Biretta, J.A., Junor, W., Perlman, E.S., Sparks, W.B., Wilson, A.S. 2003, ApJ, 586, L41 
Haswell, C.A., Tajima, T., Sakai, J.J. 1992, ApJ, 401, 495 
Junor, W., Biretta, J. A., Livio, M. 1999, Nature, 401, 891 
Junor, W. & Biretta, J.A. 1995, AJ, 109, 500 

Le Bohec, S. and the VERITAS Collaboration 2001, in Proc. 27th Int. Cosmic Ray. Conf, Hamburg, vol. 7, ed. M. Simon et 

al. (Copernicus Gesellschaft, Katlenburg-Lindau) 2643 
Le Bohec, S. and the VERITAS Collaboration 2003, in Proc. 28th Int. Cosmic Ray. Conf., Tsukuba, Japan, vol. 5, ed. T. Kajita 

et al. (Universal Academic Press, Tokyo, Japan) 2627 
Longair, M.S. 1994, High Energy Astrophysics Vol II, Cambridge University Press, Cambridge 
Lovelace, R.V.E. 1976, Nature, 262, 649 

"The MAGIC Telescope Project Technology and performance aspects": http://hegral.mppmu.mpg.de/MAGICWeb 
Marconi, A., Axon, D.J., Macchetto, F.D., Capetti, A., Soarks, W.B., Crane, P. 1997, MNRAS, 289, 21 
Melrose, D.B. 1980, Plasma Astrophysics Vol I, Gordon & Breach, New York 

Mori, M. 2000, in: "International Symposium on High Energy Gamma-Ray Astronomy", Heidelberg; 

see also http://icrhp9.icrr.u-tokyo.ac.jp/c-iii.html 
Miicke, A., Engel, R.R., Rachen, J.P., Protheroe, R.J. & Stanev, T. 2000, Comm. Phys. Comp., 124, 290 

Miicke, A. & Protheroe, R.J. 2000, in Proc. workshop "GeV-TeV Astrophysics: Toward a Major Atmospheric Cherenkov 

Telescope VI", AIP Conf Proc, Vol 515, ed. B.D. Dingus et al, 149 
Miicke, A. & Protheroe, R.J. 2001, Astropart. Phys., 15, 121 

Miicke, A., Protheroe, R.J., Engel, R., Rachen, J.P., Stanev, T. 2003, Astropart. Phys., 18, 593 
Owen, F.N., Eilek, J.A., Kassim, N.E. 2000, AJ, 543, 611 

Perlman, E.S., Sparks, W.B., Radomski, J., Packham, C, Fisher, R. S., Pina, R., Biretta, J. A. 2001, ApJ, 561, L51 



10 



A. Reimer et al.: M87 as a misaligned Synchrotron-Proton Blazar 



Prothcroc, R.J., Donca, A.-C, Rcimcr, A. 2003 , Astropart. Phys., 19, 559 

Protheroe, R.J. & Clay, R.W. 2004, Publications of the Astronomical Society of Australia, Vol. 21, 1 
Rachen, J.P. & Meszaros 1998, P., Phys. Rev. D, 58, 123005 

Rachen, J. 2000, in Proc. workshop "GeV-TeV Astrophysics: Toward a Major Atmospheric Cherenkov Telescope VI", AIP Conf. 

Proc, Vol 515, cd. B.D. Dingus et al., 41 
Reimer, O., Pohl, M., Srcckumar, P., Mattox, J.R. 2003, ApJ, 588, 155 
Reynolds, C.S., Fabian, A.C., Colotti, A. & Rees, M.J. 1996, MNRAS, 283, 873 
Reynolds, C.S., Heinz, S., Fabian, A.C., Begelman, M.C. 1999, ApJ, 521, 99 
Sparks, W. B., Biretta, J. A., Macchetto, F. 1996, ApJ, 473, 254 
Stawarz, L., Sikora, M. & Ostrowski, M. 2003, ApJ, 597, 186 
Urry, CM. & Padovani, P. 1995, PASP, 107, 803 

Weekes, T.C. and the VERITAS collaboration 2002, Astropart. Phys., 17, 221 
Wilson, A.S. & Yang, Y. 2002, ApJ, 568, 133 
Zatsepin, G.T., Kuz'min, V.A. 1966, JETP Lett., 4, 78 



Appendix A: Steady-state electron spectrum, maximum energies, synchrotron and SSC radiation 



Consider a blob which is moving relativistically with Lorcntz factor Tj and along the jet axis that is viewed from 
an observer at angle 6. In the jet frame relativistic electrons are injected into the blob of size R' in the jet frame. 
We assume that pitch angle scattering maintains a quasi-isotropic particle distribution. Our interest is to derive the 
steady-state electron spectrum. In the following all quantities are in the co-moving frame of the jet, and we omit the 
primes for simplicity. 

The number N{E)dE of particles with energy between E and E+dE in this region is governed by electron injection, 
synchrotron cooling process and particle escape on a time scale T^sc- It is described by the kinetic equation 



_d_ 

dE 

where 



EN{E) 



N{E) 



= Q{E) 



E = - 



4 (Ttc 
3 m!? 



—\ E^ = -bE^ . 
Stt J 



(A.1) 



(A.2) 



with ctt = 6.65 • 10~^^cm^ the Thomson cross section. The first term in Eq. (A.l) describes the rate of energy loss due 

to synchrotron radiation averaged over pitch-angle (because of the isotropy of the distribution) in a magnetic field B 
(in Gauss), the second term particle escape from this region at an energy independent rate = c/R'. Assuming an 
injection rate that follows a power law, Q{E) = Qo ■ E~"^, Eq. (A.l) has the solution 
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For ae = 2 and 3 the integral can be solved analytically, giving 
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(A.4) 



(A.5) 



for ae = 3, with E~^ = bTesc and Smax the maximum injected electron energy. In the case of injection powers ae 7^ 2 
or 7^ 3 we solve the integral numerically. 

The maximum energy, i^max, is limited by balancing energy gain and losses. The electrons may gain energy e.g. 
through particle acceleration. In general the acceleration time scale may be written as tacc — E/[^{E)ecB] where 
£,{E) < I may be interpreted as an acceleration rate factor. A comparison with the observations yields typically very 
low values ^(£'e,max) at the maximum energy for electrons, while in the case of protons in hadronic models the rate 
factor is much higher at the maximum proton energy, typically C(^p,max) = IO~^...l. The large difference between 
^(-Ep^max) and 5(i?e,max) Can bc naturally imderstood from the theory of plasma turbulence, since the electrons probe 
much smaller turbulence scales than the protons (Rachen (2000)). A more quantitative treatment of this issue has 
been presented by Biermann & Strittmatter (1987). 
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In acceleration theory the rate of energy gain is sensitively dependent on the upstream particle mean free path A, 
which is given by 



^^^^ = 8^d{k)k 



(A.6) 



in the small angle scattering approximation (Drury 1983) and with rg the particle's gyro-radius. The magnetic tur- 
bulence spectrum I(k) is usually expressed as a power law of the wave number k in the turbulent magnetic field: 
I{k) oc k~^. /? = 5/3 corresponds to Kolmogorov turbulence, while /? = 1 corresponds to a fully-tangled magnetic field 
resulting in "Bohm diffusion" , and is often considered for simplicity. For strong magnetic fields, Kraichnan turbulence 
f3 = 3/2 may be present. The (parallel) diffusion coefficient is then given by K|| = |A||f where v is the particle's speed 
and A|| is its mean free path parallel to the magnetic field. Hence, oc where 6 = (2-/3), and the acceleration 
time scale for the relativistic electrons and protons (rg oc E) can then be expressed by face oc oc E/^{E). In 
the following, we consider 5 to be a free parameter, and restrict our considerations to parallel shock fronts only for 
simplicity. Obviously, ^p(-E) = £,e{E) = i{E) oc E^~^ applies to both, electrons and protons. If the electron and proton 
spectra are limited by synchrotron losses, Biermann & Strittmatter (1987) found for their cutoff energy 

oc (me/m,)(^-^)/(^+^) , 1^ oc {m,/m,f^'+'^ (A.7) 



/p,inax -^p,max 

and so one expects 

oc (me/m,)^(^-^)/(^+^) . (A.8) 

S V P, max J 

The ratio of their maximum synchrotron photon energies can readily be computed to: 

(5-35)/(H-6) 



^syn,p \ 



(A.9) 



The acceleration model parameters used to calculate the SEDs of M87 can be understood for 5 ^ 0.3 which is close 
to that for a Kolmogorov turbulence spectrum. 

To obtain the synchrotron specific luminosity, for the case of no synchrotron self-absorption, as a function of 
frequency v we convolve the particle density N(E) with the synchrotron Green's function P{u, E): 

Lq{v) = j dEP{v,E)N{E) (A.IO) 



with 



V2 



V 



P(.,^) = ^eVF^-J. (A.11) 

for relativistic particles of velocity « c and an isotropically distributed magnetic field, lo^ = eB/rrie is the electron 
gyro frequency and Vc = \J\^eE'^ / {mlc') the critical frequency after pitch angle averaging. The function F{x) with 
X = u/vc can be approximated by (Melrose 1980) 



F{x) 



X I dx' K^/^ix') « 1.85a;i/3 exp(-a;) . (A.12) 



Synchrotron self-absorption will dominate the photon spectrum at low energies. The synchrotron radiation and its 
corresponding electron spectrum may be approximated by a multiple broken power law. For each part of the particle 
spectrum that is governed by a simple power law with index cte,, one can therefore use the absorption coefficient for 
synchrotron radiation of a power law electron distribution in a randomly oriented magnetic field (Longair 1994), which 
reads 

V2e3 ( ^^A"'\^B^f 3«e. + 2\ ^ / 3ae.. + 22\ ^_(„.,,+4)/2 (A.13) 



where go is Qo divided by the source volume. The synchrotron specific luminosity in the jet frame is then 

L{v) = ^ [1 - exp(-T.,,,a)] . (A. 14) 
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with Ti/,ssa = RXv,ssa- Note, that for high magnetic fields that are typical for the SPB-model, the electron spectrum, and 
thus its corresponding synchrotron spectrum above the synchrotron-self absorption break energy, is often completely 
determined by synchrotron losses for typical blazar 'blob' sizes. 

This synchrotron component represents the target photon field for photopion production and cascading in the 
SPB-model, and simultaneously manifests itself as the 'synchrotron hump' in the blazar SED after transformation of 
the luminosity L{v) into the observer frame. To save CPU-time wc fit this target photon field with a multiple broken 
power law, which is then used as an input into the SPB Monte-Carlo code (see Sect. 3). 

The importance of Inverse Compton scattering off the synchrotron photons produced by the same primary electron 
component is determined by the ratio of the synchrotron photon energy density and the magnetic field energy density. 
For B > 0.6(wphot/10^°eV cm~^)^/^ Gauss the target photon density Wphot is smaller than the magnetic field energy 
density, which is typically true for hadronic SPB models, and inverse Compton losses of the primary electron population 
can usually be neglected. In the following wc calculate the expected SSC (jet frame) specific luminosity from an electron 
spectrum Eq. (A. 3) in the Thomson limit which is valid for E <C (mgC^)^ /hv. In the ^-function approximation it is 
given by: 



where 7e,min"^eC^ is the minimum injected electron energy and n(e) is the (jet frame) synchrotron photon density. For 
the calculations in Fig. 1 and 3 we performed the e-integration numerically. 




(A.15) 



